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The propagation of an electromagnetic pulse in a plasma is studied for pulse durations that are comparable to the 
plasma period. When the carrier frequency of the incident pulse is much higher than the plasma frequency, the pulse 
propagates without distortion at its group speed. When the carrier frequency is comparable to the plasma 
frequency, the pulse is distorted and leaves behind it an electromagnetic wake. 
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I. INTRODUCTION 



II. ANALYSIS 
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The propagation of an electromagnetic wave in a 
medium Q is controlled by the dielectric function, wich 
characterizes the response of the medium to the applied 
electromagnetic held. The dielectric function of a plasma 
is 1 — ujp/uj 2 , where uj p is the plasma frequency, and ui 
is the frequency of any Fourier component of the wave. 
This simple formula also characterizes the response of 
a dielectric medium when the Fourier spectrum of the 
wave contains frequencies that are much higher than the 
resonance frequencies of the medium. 

When a monochromatic wave of frequency u> is incident 
upon a vacuum-plasma boundary, a fraction 2ki / (hr + 
kj) is transmitted and a fraction (hp — kr)/{kr + &i) is 
reflected, where ki = to / 1 c is the wave number of the inci- 

1/2 

dent wave, and fcy = (to 2 — u> 2 ) j c is the wave number 
of the transmitted wave. Now consider an electromag- 
netic pulse with carrier frequency oj c and envelope fre- 
quency uj e . The formulas for the transmission and re- 
flection of a monochromatic wave are also valid for a 
long pulse, provided one substitutes uj c for w. When 
lu c < uj p , the incident pulse is reflected completely. When 

lu c > uip, the transmitted part of a long pulse propagates 
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without distortion at its group speed c(l — oj 2 /uj 2 ) . 
Eventually, the transmitted pulse disperses. These re- 
sults are known to be valid for oj e <C w p . In this paper 
we study electromagnetic propagation in the complemen- 
tary regime oj e ~ oj p . Short-pulse propagation is gen- 
erally relevant when the long-envelope approximation is 
not valid. A specific example is the wakefield accelerator 
concept H,||. 

We use Laplace transform and Green function tech- 
niques to analyze the interaction between the laser pulse 
and the plasma. We find that the interaction can be 
divided into two stages, one in which temporal transmis- 
sion and reflection occurs at the vacuum-plasma bound- 
ary, and one in which the transmitted and reflected pulse 
propagate in the plasma and vacuum, respectively. We 
then present details of what happens at each stage, for 
incident pulses of varying carrier frequency and duration. 



We consider a laser pulse with electric field E(t, x) that 
propagates in vacuum when x < 0, enters the plasma at 
x = 0, and propagates through the plasma for x > 0. 
We assume that the plasma is characterized by a plasma 
frequency u) p . The wave equation obeyed by E(t, x) is 
given by 



c 2 d 2 xx +tu 2 p )E(t,x) = Q, 



(1) 



where df t and d 2 x 



are second-order partial derivatives 
with respect to time t and space x, and where c is the 
speed of light in vacuum. 

Let Lu p t —> t, tUpx/c — * x, so that t and x become 
dimensionless. Then the wave equation (111) becomes 



(Aft 



dl 



l)E(t,x) 



0. 



(2) 



In general, some fraction of the incoming laser pulse is 
reflected at the vacuum-plasma boundary, while the rest 
is transmitted into the plasma. We denote the incident 
electric field by Ej(t,x), the reflected field by Eji(t,x), 
and the transmitted field by Ex(t,x). Since the electric 
field is continuous across the boundary El, 



E T (t,0) + E R {t,0) =Er(t,0). 



(3) 



Similarly, since the magnetic field of the pulse is contin- 
uous across the boundary O, 



d x E I (t,0) + d x E R (t,0) = d x E T {t,0). 



(4) 



We next take the temporal Laplace transform of Eqs. 
(0) and (|i|) to obtain the equivalent boundary conditions 
in Laplace space, 

Ei(s,0) +E R {s,0) =E T (s,0), 
d x E I (s,0)+d x E R (s,0) = d x E T {s,0). (5) 

In general, the incident field Ej(t,x) propagates to 
the right (toward the plasma), while the reflected field 
E R (t,x) propagates to the left (away from the plasma). 
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We may therefore assume that Ej(t, x) and E R (t, x) have 
the space-time dependencies 

E I (t,x)=E I (t~x), 

E R (t,x)=E R (t + x), (6) 

which are consistent with the reduced equations 

{d t + d x )E I (t,x)=0, 
{d t -d x )E R {t,x)=0. (7) 

By taking the temporal Laplace transform of Eqs. @), 
and letting x — > 0, we obtain the boundary expressions 



d x E R (s,0) = sE R (s,0). 



(8) 



The Laplace transform Et(s,x) of the transmitted 
field Ex{t,x) satisfies the equation 

[dl x -(s 2 + l)]E T { S ,x) = 0, (9) 

which follows from (|J). We choose the causal solution 
(note that x > 0) 

E T (s,x) = E T (s,0)exp[-(s 2 + l) 1/2 x], (10) 

so that, at the boundary x = 0, we have 

d x E T {s,0) = -(s 2 + l)V 2 E T {s,0). (11) 

Substitution of (|J) and ( [TT|) into (||) yields the bound- 
ary condition 



sEi(s, 0) - sE R (s, 0) = (s 2 + l) 1/2 E T (s, 0). 
Equations (||) and (|l^) imply that 
— s - fs 2 + lW 2 — 



(12) 



S+^ + l) 1 



^t(s,0) 



g 2 2 l ^1/2 ^(^.0)- 



S+^ + l) 1 

It follows from the second of Eqs. (0) that 

E R (s, x) = Er(s, 0) exp(sx). 
Finally ©, ©, and (||) yield 



(13) 



(14) 



— s - (s 2 + lW 2 — 

E R (s,x) = . exp(sx)Ei(s,0), 



1)V2 



E T (s,x) 



2s 



S + (s 2 + l)V 2 

x £j(s,0). 



exp 



[-(s 2 + l) 1/2 z] 



(15) 



The coefficients of Ei(s,0) in Q15| ) are just the Green 
functions T R (s,x) and Ly(s,a:) in Laplace space for the 
reflected and transmitted pulse, respectively. We write 
the reflection Green function in the form 



T R (s,x) = R(s)G R (s,x), 



(16) 



where 



R(s) 



(s 2 + 1)V2 



S+ (S 2 + 1)1/2' 

G R (s, x) = exp(sir). 



(17) 



From the above discussion, it is clear that R(s) represents 
the reflection of the incident pulse at the vacuum-plasma 
surface, whereas the factor G R (s, x) accounts for the sub- 
sequent propagation of the reflected pulse in vacuum. 

Similarly, we write the transmission Green function in 
the form 



T T (s,x) = T(s)G T {s,x), 



(18) 



where 



T(s) = 



2s 



S+(S 2 + 1)V2' 

G T (s,x) =exp[-(s 2 + l) 1 / 2 a 



(19) 



Here T(s) represents the transmission of the incident 
pulse across the vacuum-plasma surface, whereas the fac- 
tor Gt (s, x) represents the subsequent propagation of the 
transmitted pulse in the plasma. 

We see from (0) and @ that R(s) and T(s) are 
related through the equation 



T{s) = 1 + R{s), 



(20) 



which just states the fact that the electric field is con- 
served. 

The influence of pulse duration and carrier frequency 
on the pulses' transmission and subsequent propagation 
in a plasma can be investigated by considering boundary 
fields Ei(t,0) of the form 



Ei(t,0) = exp(-w 2 t 2 )cos(w c f). 



(21) 



The parameters uj e and lu c are measures of the temporal 
envelope width and carrier frequency respectively, of the 
incident pulse at the x = boundary. We give in Table | a 
classification of the incident pulses (Ell) at the boundary. 



TABLE I. Definition of the pulse classification scheme em- 
ployed in the text. 



Pulse characteristic 



Parameter regime 



Long duration (LD) 
Intermediate duration (ID) 
Short duration (SD) 
Low frequency (LF) 
Intermediate frequency (IF) 
High frequency (HF) 



UJe < 1 
UJ e ~ 1 
LOe > 1 
LO c < 1 
LOc — 1 
Ulc > 1 
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The inverse temporal Laplace transform of (|2 
T(t) = S(t) - (2/t)J 2 {t)H(t). 



(22) 



T(t) represents the part of the laser-plasma interaction 
in which the incident pulse is transmitted across the 
vacuum-plasma boundary x = 0. The hrst term in (E^) 
represents the undistorted transmission of a pulse into 
the plasma, while the second term represents the reflec- 
tion R(t) at x = 0, 



R(t) = -(2/t)J 2 (t)H(t). 



(23) 

This is evident by comparing ( pp| ) with (gg). Equa- 
tion (|23|) shows that the reflection of the laser pulse at 
the vacuum-plasma boundary is not instantaneous, but 
rather a decaying, oscillatory function of time. This in- 
dicates that there is a harmonic response in the plasma 
to the incident pulse, which produces a delayed, rather 
than instantaneous, reflected pulse. This response takes 
the form of harmonic oscillations of plasma charges about 
their equilibrium positions, which are induced by the in- 
cident sinusoidal pulse. 

One can investigate the dependence of the reflected 
pulse at x = on the duration and frequency of an im- 
pinging pulse Ei(t,0) by calculating the convolution 



E R (t,0) 



Ei(t', 0)R{t-t')dt' 



(24) 



for different values of the parameters uo e and lo c in 
Ei(t,0), where £r(*,0) is given by @. Figure [j] shows 
the reflection response for incident pulses of intermedi- 
ate duration (ID), with carrier frequencies varying from 
intermediate (IF) to high (HF). Figure || shows the reflec- 
tion response for incident pulses of short duration (SD), 
again with carrier frequencies varying from intermediate 
(IF) to high (HF). It is seen in Figs. [I] and g that the 
reflection response diminishes as the carrier frequency of 
the pulse is increased. We also note that as the duration 
of a pulse is shortened (i.e., as Lo e is increased beyond 
1), the reflection response diminishes. This is consistent 
with the fact that, as an incident pulse is shortened, more 
of it will already have entered and propagated into the 
plasma before the plasma's delayed reflection response 
[as described below (p3|)] takes place. In particular, if 
u> e 3> 1, the pulse is transmitted completely, with no 
distortion. 
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FIG. 1. Temporal evolution of the reflection response 
En(t,0) [Eq. (E4)] at the vacuum-plasma boundary, for an 
incident pulse of intermediate duration (u e = 1) and carrier 
frequency ui c = 1 (dotted), lo c — 3 (dashed), and uj c — 10 
(solid). 
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FIG. 2. Same as in FigjlJ but for an incident pulse of short 
duration (uj e — 5). 

The propagation of the reflected pulse in vacuum is 
characterized by the function Gfi(t, x), which is the in- 
verse of Gr(s 7 x) in (|i~7|), 



G R (t,x) = 5(t + x). 



(25) 



This means that the reflected pulse Eji(t,x) has the 
space-time dependence Eji(t,x) — Eji(t + x), and prop- 
agates through the vacuum in the negative s-direction 
away from the vacuum-plasma boundary, and without 
distortion. 

We now focus on the transmitted pulse ET(t,x). The 
propagation of the transmitted pulse through the plasma 
is characterized by the function GT(t,x) given by the 
inverse of Gt(s,x) in ([l9|). Gr{t,x) is found by first 
writing Gt(s,x) as the spatial derivative 



where 



G T (s,x) = -d x F T (s,x), 



( s 2 + 1 )l/2 

The inverse of Ft(s, x) is given by Q| 

F T (t,x) = J [(t 2 -x 2 ) 1/2 }H{t-x), 

so that 

G T {t,x) =S(t-x) 

■h[{t 2 -x 2 y/ 2 ] 

X — 777~T77. H (t — X) 



(26) 



(27) 



(28) 



(t 2 



r 2)l/2 



(29) 



Equation ([29|) represents the combined effect of a distor- 
tionless propagation of the transmitted pulse (first term) 
and the propagation of a dispersive wake generated by 
the plasma (second term). 
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We next compute the total Green function x) by 
inverting (|l8|). One way to do this is to compute Tx{t 7 x) 
as the convolution 



r T (t,x) 



T{t-t')G T (t',x)dt', 



(30) 



where T(t) is given by @, and G T {t, x) by @. Again, 
( |30| ) clearly shows the two-stage process of transmission 
followed by propagation. Analytic evaluation of (^0|) is 
quite involved. However, there is a simpler method for 
obtaining Tx{t,x) analytically that avoids integration, 
and requires only the computation of derivatives. From 
( |l8| ) and (|l9|), we see that IV (s, x) can be written as the 
derivative 



where 



T T (s,x) = -2d x [sf(s,x)], 
exp[-(s 2 + l) 1 /^] 

[s+(.5 2 + l) 1 /2]( s 2 + 1)1/2 

f(s,x) has the inverse transform Q 

' t 



f(s,x 



(31) 



(32) 



1/2 



F(t,x) = T(t,x)H(t - x) 



x 



M(t 2 - 



t + X 

2 )V 2 ]ff (t - a:). 



(33) 



This inverse transform only holds for x > 0, which is in 
accord with our assumptions of the pulse entering the 
plasma at x = 0, and propagating into the plasma for 
x > 0. Since F(0 + ,a;) = for x > 0, we have from 
standard Laplace transform theory that dtF(t, x) is the 
inverse transform of sf(s,x). Therefore, from (pl|), we 
have 



r T (i,x) = -2d t y^,z) 



(34) 



The term —2d 2 x J-(t, x) in ([34]) represents a modification 
to the incident pulse, caused by reflection at the vacuum- 
plasma boundary and dispersion in the plasma. It is 
given by 



2dlT(t,x) = - 




,t + x ' (t + x) 2 
x J 2 (t 2 -x 2 ). 



(35) 



From (p"5|), JIg|), and (|l9|), we see that the transmitted 
pulse Exit, x) is given by the Green function integral 



E T (t,x) 



Ei(t',0)T T (t-t',x)dt'. 



(36) 



We next perform the integration in ( pq ) for incident 
pulses of the form (21). We first consider incident pulses 



of intermediate duration. Figures || and || show plots 
of the propagation of ID-HF and ID-IF incident pulses, 
respectively. 
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FIG. 3. Spatial dependence of a transmitted pulse 
E T (2Q,x) [Eq. @] at time t = 20 (solid). The incident 
pulse crossed the vacuum-plasma interface at t — 0, and had 
a spatial dependence in vacuum characterized by Eq. (pl[), 
oj e = 1 (ID), and u c = 10 (HF). The incident pulse's spatial 
dependence translated to t — 20 is shown by the dotted curve, 
for comparison with the resulting transmitted pulse. 
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FIG. 4. Same as in FigJ^, but with incident pulse parame- 
ters w e = 1 (ID) and u c = 1.5 (IF). 

We see that the high frequency (HF) incident pulse 
propagates practically undisturbed across the vacuum- 
plasma interface and into the plasma, while the 
intermediate-frequency pulse develops an electromag- 
netic (EM) wake. In the Appendix, we derive the fol- 
lowing perturbative expansions for v g and v p in the high- 
frequency case: 



2 /2 - e 4 /i 

2 /o i o,4 , 



1/vg Pa 1 + e 2 /2 + 3e 



(37) 



where e = lo p /lo c . The right sides of Eqs. ( p7|) are just 
the first three terms in the MacLaurin expansions of (1 — 
e 2 ) 1 ' 2 and (1 - e 2 )- 1 ' 2 . 

We next consider the propagation of short (SD) inci- 
dent pulses. Figure || shows a plot of an incident SD-IF 
pulse. 



4 



LU 



0.06 
0.04 
0.02 


0.02 
0.04 
0.06 




16 17 18 19 20 21 

x 



FIG. 5. Same as in Fig^, but with incident pulse parame- 
ters Lo e = 5 (SD) and lo c = 1.5 (IF). 

As expected, the wake generation is smaller than for 
the incident ID-IF pulse. And as the frequency of the 
incident SD pulse is increased, it is found that wake gen- 
eration is practically nonexistent. 



III. SUMMARY 

In this paper, we considered the transmission and re- 
flection of an electromagnetic pulse at a vacuum-plasma 
boundary, and the subsequent propagation of the trans- 
mitted pulse in the plasma. We extended the well-known 
theory for long pulses into the short-pulse regime, in 
which the pulse duration is comparable to the inverse 
plasma frequency. When the carrier frequency of the in- 
cident pulse is much higher than the plasma frequency, 
most of the incident pulse is transmitted without dis- 
tortion. Subsequently, the transmitted pulse propagates 
without distortion at its group speed. When the car- 
rier frequency is comparable to the plasma frequency, 
the transmitted pulse is distorted, and leaves behind it 
an electromagnetic wake. The reflected pulse is delayed 
relative to the incident pulse, and is also distorted. When 
the carrier frequency is less than the plasma frequency, 
the incident pulse is absorbed by the plasma before being 
reemitted. 
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APPENDIX A: PROPAGATION OF A 
HIGH-FREQUENCY PULSE 

Let bj c t — > t, uj c x/c — > x, and w p /w c — > e, so that t and 
x become dimensionless. Then the wave equation (Q) can 
be written as 



(4 - dl 



e )E = 0. 



(Al) 



The study of pulse propagation is facilitated by the char- 
acteristic transformation 



= t — (3x, £ = x — /3t, 



(A2) 



where (3 < 1. In terms of the characteristic variables r 
and £, the wave equation (Al) can be rewritten as 



[{l-f)(dl T -dl) + ^]E = Q. 



(A3) 



One can solve Eq. (A3) by using multiple scale analysis 
g. To do this, one introduces the time and distance 
scales 



r„ = e"r, = e»f. 
Correct to second order, one can write 



Or 

9( 



d Tn + ed T , + t 2 d T 



% + e% 



e 2 %. 



(A4) 



(A5) 



Guided by the well-known characteristics of a long pulse, 
we assume that 







1 



e 2 /3 2 



and 



E(r,0 = B(r 2 ,a)exp(-iro). 



(A6) 



(A7) 



Ansatz ( A7) corresponds to a pulse that has a carrier fre- 
quency of unity and an amplitude that varies on the slow 
scale £i. For this amplitude variation, /3 is the group 
speed of the pulse, and the characteristic variables are 
proportional to time and distance me asu red in t he pulse 
frame. One now substitutes Eqs. (A5) - (A7) in Eq. 
( [A3] ) and collects terms of like order. The zeroth- and 
first-order equations are satisfied automatically by con- 
struction. 

In second order, 



It follows from Eq. (A8) and ansatz (A7) that 
02 = -1/2. 

In third order, 

-4p 2 (d*-dl £l )E = 0. 



(A8) 



(A9) 



(A10) 
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Equation (AlC) is consistent with ansatz (A7), in which 
E is assumed to be independent of £o and t± . In fourth 
order, 



[4^ oT2 -9| o6 ) + 2/3 2 (^ lTl -a| l6 ) 



9 

T T 2 

+ (2/3 4 + fi)(d 2 T0T0 - a| 0?0 )] £ = 0. (All) 



The pulse has a carrier frequency of unity by construc- 
tion, so the dependence of E on r 2 cannot be oscillatory. 
It follows from this constraint that (2/34 + /?!) = and, 
hence, that 



ft = -1/ 



(A12) 



The group speed fi ~ 1 — e 2 /2 — e 4 /8, which is just the 
first three terms in the Maclaurin expa nsion of (1 — e 2 ) 1 / 2 . 
The remaining nonzero terms in Eq. ( All ) are 



(2id T2 +dl A )B = Q, 



(A13) 



which describe the dispersal of the pulse. 

Finally, note that ansatz (A7) constrains the phase 
speed to be the inverse of the group speed. Since no 
contradictions appear in the s ubse quent analysis, the as- 
sumptions underlying ansatz (A7) are correct. One can 
also use the ansatz 



E(t,0 = B(r 2 ,£i)exp[^o - i(l - ^)r ] 



(A14) 



which does not constrain the phase speed, but leads to 
the same result. 
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